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SOME  CONICAL  AND  QUASI -CONICAL  FLOWS  IN  LINEARIZED 
SUPERSONIC -WING  THEORY 
By  Herbert  S .  Ribner 


SUMMARY 

A  number  of  conical  and.  quasi -conical  linearized  supersonic 
flows  have  been  derived „  These  flows  may  be  applied  in  lift- 
cancellation  techniques  in  the  determination  of  wing-lift  dis¬ 
turbances  that  arise  at  subsonic  trailing  edges.  Some  of  the 
results  are  applied  to  damping  in  roll  and  pitch  in  another  paper. 

Two  methods  of  analysis  have  been  employed,  both  involving 
integral  equations.  One  is  a  development  of  the  membrane  method 
of  Eward;  the  other  is  an  analog  of  the  superposition  method  of 
Schlichting.  In  both  cases,  source  distributions  are  used.  The 
mathematical  equivalence  of  the  first  method  and  a  new  doublet - 
distribution  method  of  Goodman  and  Mirels  is  shown. 


INTRODUCTION 

The  analysis  of  the  flow  over  a  sweptback  wing  (reference  1) 
makes  use  of  Lagerstrom's  concept  of  lift  cancellation  (reference  2). 
The  starting  point  is  a  delta  wing  of  infinite  chord.  Suitable  flows 
are  superimposed  that  cancel  the  lift  outboard  of  and  behind  a 
certain  boundary.  This  boundary  is  chosen  to  constitute  the  tips 
and  the  trailing  edge  of  the  sweptback  wing. 

The  cancellation  flows  for  the  wing  tip  modify  the  lift  in  the 
tip  region;  they  are  not  of  concern  herein.  The  cancellation  flows 
for  the  trailing  edge  modify  the  lift  in  a  region  ahead  of  the  trail¬ 
ing  edge  if  the  edge  is  "subsonic "(that  is,  if  the  component  stream 
velocity  normal  to  the  edge  is  subsonic).  Now  the  lift  distribution 
of  a  delta  wing  at  an  angle  of  attack  is  substantially  flat  in  the 
center.  (See  fig.  1.)  Thus,  an  approximate  cancellation  behind  the 
trailing  edge  would  be  afforded  by  a  constant  lift  such  as  flow  I  in 
the  figure.  (This  idea  was  originated  in  reference  1.)  .Correspond¬ 
ing  approximate  cancellation  flows  (ill  and  TV)  for  rolling  and 
pitching  motions  are  also  shown  in  the  figure.  The  derivation  of 
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these  three  flows  is  the  main  object  of  this  report;  the  work  was  > 

initiated  at  the  MCA  Langley  laboratory  and  completed  at  the  NACA 
Lewis  laboratory* 

The  analysis  leads  to  integral  equations,  some  with  several 
apparent  solutions .  The  choice  of  the  physically  correct  solution 
in  each  case  is  determined  by  comparison  with  the  solution  for  a 
closely  related  flow  that  is  more  easily  calculated.  Thus, 
altogether,  derivations  are  given  for  eight  flow  distributions. 

Four  of  these  flows  are  conical:  the  velocity  components  and  pres¬ 
sure  are  constant  along  any  ray  from  the  vertex.  The  remaining 
four  are  quasi -conical:  the  velocity  components  and  pressure  are 
proportional  to  the  distance  along  any  ray. 

The  practical  application  of  several  of  these  cancellation 
flows  in  the  evaluation  of  damping  in  roll  and  pitch  is  made  in 
reference  3. 


CANCELLATION  FLOWS 

Each  of  the  eight  flows  discussed  can  be  considered  to  repre¬ 
sent  a  slightly  cambered  lifting  surface  lying  essentially  in  the 
z  =  0  plane.  (All  symbols  are  defined  in  appendix  A.)  Thus,  the 
u  and  v  velocities  are  antisymmetric  with  respect  to  the  z  «  0  plane, 
and  the  w  velocity  is  symmetric.  With  these  reservations,  the 
boundary  conditions  can  be  specified  for  the  upper  surface  only. 

(This  simplification  is  employed  throughout  the  report.) 

The  flows  contemplated  for  the  partial,  cancellation  of  the 
lift  behind  the  trailing  edge  are  schematically  designated  in 
figure  2  for  the  several  cases .  A  particular  cancellation  flow 
(any  one  of  I  to  IV,  fig.  2)  is  superposed  on  the  sweptback  wing 
so  that  the  w  =  0  regions  lie  ahead  of  the  trailing  edge  and  the 
shaded  region  lies  behind.  The  lift,  and  hence  the  u  velocity,  is 
specified  in  the  shaded  region,  and  this  specification  determines 
the  lift  cancellation.  (See  fig.  1.)  The  specification  w  =  0  for 
the  part  -of  the  flow  overlapping  the  wing  ensures  that  the 
resultant  flow  shall  be  tangent  to  the  surface  of  the  wing;  the 
basic  delta-wing  flow  already  provides  the  correct  value  of  upwash  w 
on  the  wing  (for  example,  -<xV  for  angle  of  attack)  and  the  cancel¬ 
lation  flow  must  therefore  add  none . 

The  u  velocity  of  the  cancellation  flow  in  the  regions  w  =  0 
is  algebraically  additive  to  that  of  the  basic  delta  flow  there. 

This  additive  u  velocity  corresponds  to  a  pressure  disturbance  caused 
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by  the  wing  trailing  edge.  The  central  problem  of  this  report  is 
the  determination  of  the  u  velocity  for  each  of  the  cancellation 
flows  I  to  IY  from  the  boundary  conditions  stated  in  the  preced¬ 
ing  paragraph. 

In  the  integral -equation  method  of  solution  as  applied  herein, 
several  solutions  may  be  found  for  certain  of  the  integral  equ¬ 
ations.  Unique  solutions  can  be  obtained  in  each  case,  however, 
for  a  flow  that  is  closely  related  to  the  desired  flow.  The 
correct  solution  for  the  desired  flow  may  be  identified  by  a  com¬ 
parison  of  the  nature  and  the  location  of  singularities  in  the  two 
flows.  The  identification  is  made  more  convincing  by  a  quanti¬ 
tative  comparison  of  the  u  velocity  over  the  region  of  interest  in 
addition  to  a  consideration  of  the  singularities. 

Thus  in  figure  2  flows  I  to  IV  are  the  desired  flows  and  flows 
I'  to  IV '  are  the  respective  related  flows.  The  boundary  conditions 
for  the  desired  and  related  flows  differ  only  in  the  presence  or 
absence  of  a  left-hand  w  =  0  region.  Because  of  the  relative  remote¬ 
ness  of  the  left-hand  region,  conditions  in  this  region  may  be 
expected  to  affect  but  slightly  the  u  velocity  induced  in  the  right- 
hand  w  =  0  region,  so  long  as  that  region  is  relatively  narrow.  With 
this  limitation  the  u  velocity  in  the  right-hand  w  =  0  region  cal¬ 
culated  for  flow  I ' ,  then,  is  presumed  to  be  a  good  approximation  to 
that  for  flow  I,  and  similarly  for  the  other  pairs  of  flows. 

It  will  be  convenient  to  obtain  first  the  solutions  for  the 
related  flows  I'  to  IV'  so  that  they  will  be  available  to  aid  in 
identifying  the  correct  solutions  for  flows  I  to  IV,  respectively. 
Flows  I'  to  IV'  require  only  a  simple  inversion  of  an  Abel-type 
integral  equation.  Flows  I  to  IV  require  a  more  elaborate  procedure. 


Flow  I ' 

(For  Angle  of  Attack) 

In  flow  I'  (fig.  3)  the  w  velocity  (or  surface  slope)  is 
prescribed  in  region  A,  but  not  the  u  velocity  (or  lift).  The 
u  velocity  is  prescribed  in  region  B,  but  not  the  w  velocity.  If 
the  w  velocity  were  known  in  region  B,  this  wing  would  be  of  pre¬ 
scribed  camber.  According  to  Puckett  (reference  4),  the  wing 
could  be  represented  by  a  distribution  of  sources  in  proportion  to 
the  local  value  of  w.  The  unknown  u  velocity  in  region  A  could 
then  be  obtained  by  simple  integration  and  differentiation.  The 
source  representation  can  still  be  used,  however,  even  though  the 
distribution  of  w  in  region  B  is  unknown.  In  this  case  the  pre¬ 
scribed  conditon  on  u  in  region  B  will  give  rise  to  a  soluble 
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integral  equation  for  the  unknown  w  distribution.  This  procedure 
is  a  development  of  Evvard's  method  (reference  5).  The  details 
follow. 


The  surface  u  velocity  is  given  in  reference  5,  equation  (3a). 
In  slightly  modified  form,  and  in  the  present  notation,  it  is 


u(£,n)  =  - 


(i) 


The  oblique  coordinates  of  reference  5  (Mach  coordinates)  are  employed. 
(See  fig.  3.)  Unsubscripted  values  refer  to  the  field  point,  sub¬ 
scripted  values  to  the  source  point.  The  transformations  relating 
(x,y)  and  (|,tj)  are 


x  =  m 
M 


t  -  §  <*-<*> 


By  =  §  Cn-S)  n  »  jfJ-  U+By) 


(2) 


The  surface  integral  is  taken  over  the  area  S  bounded  by 
the  forward  Mach  cone  frcan  the  point  (|,tj).  The  line  integral  is 
taken  along  any  lines  within  S  or  bounding  S  across  which  the 
w  velocity  experiences  a  jump  Aw.  The  |  -axis  and  the  line  q  =  T  { 
(fig.  3)  might  be  such  lines.  Consider,  however,  the  requirement 
of  general  similarity  with  flow  I.  The  imposition  of  the  Kutta 
condition  at  the  trailing  edge  of  the  wing  in  figure  1  requires  the 
continuity  of  w  across  the  boundaries  q  =  T  £  and  |  **  TH  of  the 
shaded  region  of  flow  I,  as  sketched  in  figure  2 .  For  flow  I ' ,  w 
must  similarly  be  continuous  along  the  right-hand  edge  of  the  shaded 
region  (t)  =  t£).  Also,  the  continuity  of  w  across  the  £-axis 
(Mach  line)  can  be  shown  to  be  a  consequence  of  the  finite  lift  along 
that  line.  The  line  integral  in  equation  (1)  accordingly  vanishes 
for  flow  I ' . 


Let  the  point  (|,tj)  be  located  at  P  (fig.  3);  P  repre¬ 
sents  an  arbitary  point  in  region  B.  Then  the  forward  Mach  cone 
includes  the  shaded  area  in  region  B  and  an  adjacent  unshaded  area 
in  region  A.  The  quantity  Sw/^  =  0  in  region  A;  therefore, 

only  the  shaded  area  contributes  to  the  surface  integral. 

Equation  (1)  may  thus  be  written 
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u(|,T))  =  - 


_1_ 


region  B 


(3) 


The  unknown  inner  integral,  multiplied  "by  l/«M,  may  "be  designated 
which  gives 


c  n 


region  B 


(4) 


Equation  (4)  is  an  Abel-type  integral  equation  for  the  unknown 
function  J(|,tvj_).  Its  solution  is  given  (reference  6)  by 

-Til 

aiLiliii  (5) 

0  ^ 

(The  more  general  treatment  of  Abel's  integral  equation  in  refer¬ 
ence  7  implies  that  equation  (5)  will  provide  singular  as  well  as 
continuous  solutions.  The  only  important  restriction  is  that 
u(|,T])  must  be  such  that  the  integral  on  the  right-hand  side  of 
equation  (5)  is  continuous.  In  reference  6,  however,  certain 
additional  restrictions  are  placed  on  u(t,T])  to  limit  equation  (5) 
to  continuous  solutions.  In  the  present  report,  these  additional 
restrictions  are  disregarded  and  singular  solutions  may  be  expected, 
as  in  equation  (6).) 

In  the  present  case  the  function  u(!,t])  is  a  constant  uQ 
(fig.  3).  The  solution  for  J  is  therefore 


J(Mi)  =  -7 


1  3 
* 


J(Mi) 


Now  let  the  general  point  (£,q)  be  located  in  region  A 
at  Q  (fig.  3).  The  only  change  in  the  expression  for  u 


(6) 
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(equation  (3))  is  in  the  upper  limit  for  r^.  This  upper  limit  is 

now  at  the  intersection  of  the  line  PQ  with  the  line  q  =  t|  . 
Accordingly, 

* 

u(|,ri)  =  -  - - —  1  ■  region  A  (7) 


The  value  of  J  already  obtained  in  equation  (6)  applies  here  as 
well  as  in  equation  (4) .  Upon  making  the  substitution, 


u0 
=  T 


t  1 


dq 


0 


2u 


0  .  -1 

-  sm 


region  A 


(8) 


It  will  be  convenient  to  reexpress  the  result  of  equation  (8) 
in  terms  of  the  conical  coordinate  C  =  Py/x.  Along  the  line 
r|  =  ,  the  coordinate  a  is  given  the  value  n.  The  transformation 


then  yields 


u 


(o)  = 


s  = 

h  x(1-0) ^ 

n  = 

h x(1+a) 

> 

nr 

l+n 

T  = 

1-n 

J 

2u0 

it 

sin-^ 

(l+n) (l-o) 
\j  (l-n)(l+d) 

(n  <  0  < 

(9) 


(10) 


Equation  (10)  gives  the  hitherto  unknown  u  velocity  in  region  A 
of  figure  3.  The  w  velocity  in  region  B  is  still  unknown.  This 
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w  velocity  will  not  be  needed  in  the  present  investigation,  but 
the  method  of  solution  is  given  in  appendix  B  as  a  matter  of 
interest . 


Flow  II ' 

(Special) 

The  development  leading  to  equation  (5)  in  the  treatment  of 
flow  X1  leaves  the  boundary  condition  on  u  as  yet  unspecified. 
The  development  and  the  equation  are  thus  sufficiently  general  to 
be  applied  to  any  of  the  flows  I'  to  XV’ . 

The  present  case  is  complicated  by  discontinuity  in  the 
u  velocity  across  the  x-axis  (fig.  4).  The  point  P  =  (£,r|)  may 
lie  in  either  of  the  regions  B^  or  B2 •  The  corresponding 
solutions  for  J(£>t|i)  from  equation  (5)  will  be  different;  call 
them  Jj(t,Th)  and  JgdjTi)*  respectively.  The  conditions  on 
u(£,t])  are 

u(|,ti)  =  -u0  0  <  r\  <  i  region  Bx  * 

u(l,q)  =  uQ  |  <  t)  <  T|  region  B2 


Then 
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The  results  are 


(I  >  T)l )  * 


U, 


'fi! 


T  ft  s  U0  /  1 

J2(^T)l)  =  — 


(ii) 


j 

Equation  (7)  is,  like  equation  (5),  applicable  to  any  of  the 
flows  I'  to  IV'.  For  flow  II'  the  appropriate  function 
is  in  the  range  0  <  r)  <  |  and  JgU,^)  in  the 

range  |  <  q  <  t£.  Therefore 

ni  mt 


u(!,q)  =  - 


q-q. 


r>  t| 

X) 

dqi 

It 

0  *  \J 

/\^(Ti-Tii)(Tii-U 

2u, 


it 


n  -  cos 


-i  (gr-i) N . 


FI 


region  A 


In  terms  of  the  conical  coordinate  O  =  Py/x  this  equation  is, 
after  the  transformation  equations  (9)  are  applied. 


cos 


-1 


2n-o(l+n) 

c(i Sr 


(12) 
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Equation  (12)  gives  the  u  velocity  in  the  w  =  0  sector 

(n  <  O  <  1)  of  flow  II*  (fig.  1),  which  is  zone  A  of  figure  4. 


Flow  III' 

The  appropriate  value  of  u  to  he  substituted  in  equation  (5) 
is,  from  figure  2, 


u  =  Ky 


or 


u  =  I  M) 


so  that 


\  K  S 

J(^Tii)  =  -  jtM 


(t-l)dri 


K  ^^1~  \ 

*M  'sfni 


(13) 


Insertion  of  this  value  into  equation  (7.)  gives 

ua,T!)=^  -  — —  (rj  >T|) 


^1(1  "Til) 


(14) 


which  integrates  to 


u(£,t|)  = 


2K 

itM 


(t£)  sin-1^ 


-  Jt£(tt^£) 


10 


NACA  TN  2147 


In  terms  of  the  conical  coordinate  a=  py/x  this  equation  is, 
after  the  transformation  equations  (9)  are  applied. 


u(o,x)  = 


KX 

JtP 


2 a  sin 


Vz(l+n)  (l-g)  (g-n) 
1-n 


(15) 


Equation  (15)  gives  the  u  velocity  in  the  v  =  0  sector 
(n  <  O  <  1)  of  flow  III'  of  figure  2. 


Flow  TV' 

Refer  to  figure  2  for  the  value  of  u  to  he  substituted 
in  equation  (5) : 


Thus 


u  =  Kx 
- 1 


JU,^) 


pm 


KP  3 
itM  5tJi 


/\fnpn 


«p  2t)i+£ 


Substitution  of  this  value  in  equation  (7)  gives 

A 


u(£,ti) 


Kg. 

itM 


(^1+lLdJ]i  (n  >  re) 
a/^T^T 


The  integration  yields 
n  \  2*P 

U(M)  =  Tm 


M)  '  \|T*(T1"Te) 


(16) 


^  L2.  L 
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In  terms  of  the  conical  coordinate  a  =  Py/x  this  equation  is, 
after  the  transformation  equations  (9)  are  applied. 


,  \  Kx 

u(o,x)  =  — 


This  equation  gives  the  u  velocity  in  the  w  =  0  sector 
(n  <  a  <  1)  of  flow  IV'  of  figure  1. 


Flow  I 

(For  Angle  of  Attack) 

Flows  I  to  IV  differ  from  flows  I'  to  IV',  respectively,  in 
the  specification  of  a  left-hand  sector  w  =  0  to  match  the  right- 
hand  sector  (fig.  2).  The  solution  can  no  longer  he  obtained  by 
the  simple  inversion  of  an  Abel-type  integral  equation.  Kesort  is 
therefore  made  to  a  different  integral -equation  formulation  in 
which  advantage  is  taken  of  the  conical  nature  of  the  flow .  The 
method  is  an  analog  of  the  superposition  method  of  Schlichting 
(reference  8) . 

According  to  the  considerations  developed  in  the  discussion  of 
flow  I',  flow  I  may  be  represented  by  a  suitable  source  distribution. 
At  every  point  the  required  source  strength  is  proportional  to  the 
local  value  of  w.  Thus  the  source  strength  is  zero  in  the  two 
outside  sectors.  (See  fig.  2.)  In  the  central  sector  w  is 
unknown,  but  there  is  a  condition  on  the  u  velocity.  The  source 
distribution  there  must  be  so  chosen  that  this  condition  on  the 
induced  u  velocity  is  met.  The  formulation  of  this  condition. gives 
rise  to  an  integral  equation  for  the  unknown  source  distribution. 

A  convenient  form  of  integral  equation  is  obtained  as  follows.  The 
source  distribution  over  the  region  u  =  uq  is  considered  to  be 

obtained  by  superposing  uniform  sectors  of  infinitesimal  strength 
and  different  lateral  extent  ±0-^.  Section  A-A  of  figure  5  illus¬ 
trates  such  a  distribution.  The  strength  of  an  elementary  sector  is 
given  by 

(o1  >  0) 


^  4ai 


2  sin’ 


(1+n) 

1  W  2(l+n)  (l-a)  (a-n) 

/( 1-n) 

i(i+al 

1  1-n 

(17) 
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Each  elementary  source -sheet  sector  induces  a  certain  infinitesi¬ 
mal  u  velocity  and  the  total  u  velocity  is  given  "by  the  integral. 
This  integral  must  he  equated  to  the  prescribed  u  velocity  shovn 
in  figure  5.  The  result  is  the  desired  integral  equation  for  the 
unknown  distribution  strength  3w/3o^. 

The  infinitesimal  u  velocity  induced  by  an  elementary  source - 
sheet  sector  extending  from  -O^  to  can  be  obtained  from 

reference  4  (equation  (31)).  In  the  present  notation  it  is 


Au(a)  = 


- °L_  ( 3osh-1  +  cosh-1 

|a_aL| 


(18) 


where  an  absolute  value  sign  has  been  added  to  the  denominators  of 

the  cosh"**-  terms.  With  the  absolute  value  sign,  equation  (18) 
applies  both  for  \c\  <  and  for  |a|  >  Gq;  whereas  in  reference  4 

the  two  cases  are  separately  covered  in  equations  (31)  and  (33), 
respectively.  (The  same  expression  results  from  the  addition  of  two 
oppositely  swept  line  sources  of  the  acceleration  potential.  See 
reference  9,  equation  (12).) 

The  total  u  velocity  induced  by  the  superposed  source -sheet 
sectors  is 


U  (a) 


_1_ 

itP 


r  n 


cHr  al 

\Ji-g-l2 


da-j^ 


(19) 


If  u(O)  is  put  equal  to  the  specified  value  Uq  in  the  range 
| o |  <  n,  then  equation  (19)  represents  an  integral  equation  for 
the  derivative  5w(o1)/9ai  of  the  unknown  source  strength.  This 


equation  as  it  stands  is  too  complex  to  be  useful.  Great  simpli¬ 
fication  results,  however,  upon  differentiating  both  sides  with 
respect  to  o: 


r'Ti 


3u 

3a 


2a 


rtP  V 


'l-o 


35^  aid0i 

o*-<? 


(20) 


s 
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f 


♦ 


(The  Cauchy  principal  value  is  to  he  taken  for  the  integral.)  In 
the  range  |o|  <  n,  5u/Sa  =  0.  The  following  less  formidable 
integral  equation  Is  thus  obtained: 


n 


dw 

3o^Ql 

a^-o2 


cL^  =  o  |a|  <  ri 


(21) 


2  2 

The  substitution  t^  =  ,  t  =  0  reduces  equation  (21)  to  a 

special  case  of 


4' 


ni> 


'-'a 


t2~t 


=  g(t) 


(21a) 


The  integral  equation  (21a)  is  well  known  from  incompressible  thin- 
airfoil -section  theory,  but  the  usual  inversions  yield  only  the 
trivial  solution  f (t2)  =  0  for  g(t)  =0.  A  nontrivial  solution 

of  equation  (21)  for  cHr/dc^  has  been  suggested  by  C .  E.  Brown  of 
the  NACA  Langley  laboratory.  Put 


o2  =  (1  -  cos  6) 

9  r,2 

0*2  =  \  (!  -  cos  e) 


}  (22) 


=  Cf(S)  (C  =  unknown  constant) J 
Then  equation  (21)  becomes 


pjt 

f(Q)  sin  e  d8  =  0  (23) 

cos  0  -  cos  5 

-  0 


% 


* 
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Now  according  to  reference  10 

n* 


ae 


C03  9  -  cos  & 


=  0 


whence  a  solution  of  equation  (23)  is 

f(0)  =  esc  9 

dw 


A 


or 


_ —  =  C  esc  9 
oo^  j 


(24) 


The  question  of  uniqueness  is  deferred  until  later;  in  the  meantime 
equation  (24)  is  considered  as  the  physically  correct  solution. 


Equation  (24)  carries  the  solution  for  dw/do^  far  enough  to 

enable  the  determination  of  u.  Thus,  equation  (20)  is  still  valid 
in  the  range  n  <  |al  <_1.  In  this  range  du/da  ^  C,  and 
equation  (20)  may  serve  for  the  evaluation  of  u.  The  transformation 
from  Oi  to  9  (equations  (22))  is  again  convenient.  (The  trans¬ 
formation  from  a  to  5  is  not  used  as  it  leads  to  cos  6  >  1.) 

The  result  is 


du  =  CC 

55 


Jo 


9}  Bin  9  dO 
-  -  1  +  cos  9 


n  <\a\  <  1 


and  with  f (©)  =  esc  0 


Su _  Co  . 

30 


r 

d0 

Tz  " 

20  ,  -  1  +  COS  0 


n<|a|<  1 


(25) 


0 


4. 


L  V£  L 
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* 


11  (table  64,  equation  (12) ) . 


n  <  |cr|  <  1  (26) 


For  purposes  of  integration,  O  may  be  restricted  to  the  positive 
range  so  that  a  /  |a  |  =1;  the  resultant  integral  will  be  applied 
for  both  the  positive  and  the  negative  value  of  O  because  of  the 
symmetry  of  the  flow  in  O.  The  elliptic-function  substitutions  of 
appendix  C  are  helpful  in  this  and  other  more  difficult  integrations 
of  the  same  nature  in  this  report.  The  integral  of  equation  (26) 
between  the  limits  O  and  1  is  found  to  be 

u(0)  =  g-jy-  F(cP,k) 

where  F  is  the  incomplete  elliptic  integral  of  the  first  kind 
with  modulus 


This  integral  is  evaluated  in  reference 
The  result  is 


3u 

3a 


Cn2a 


2P  |a|  Jtl-o2) (a2-n2) 


and  amplitude 


cp 


The  constant  C  is  established  by  the  condition  that  u(n)  =  uQ 
so  that  finally 


u(o) 


F(cp,k) 


(27) 


Equation  (27)  is  a  solution  for  the  u  velocity  in  the  out¬ 
board  sectors  n  <  O  <  1  of  flow  I  in  figures  2  and  5.  Equation  (27) 
is  not  a  unique  solution  of  the  integral  equation  (21),  and  it  remains 
to  be  shown  that  it  is  the  physically  correct  solution  for  flow I. 

The  singularities  in  du/da,  equation  (26),  may  be  compared  with 
those  in  the  corresponding  result  for  flow  I',  obtained  by  differen 
ating  equation  (10) .  In  both  cases  half -order  singularities  are 
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found  at  o  =  n+  and  a  =  1  (that  is,  |du/do|  approaches  infinity 

as  (O’ -n) /2  when  O  approaches  n  from  above  (designated  by 
n+)  and  there  is  a  similar  behavior  as  C  approaches  1).  The 
quantitative  comparison  of  equation  (27)  with  the  corresponding 
solution,  equation  (10),  for  the  right-hand  sector  of  flow  I' 

(figs.  2  and  3)  is  even  more  convincing.  The  two  solutions  are 
plotted  together  in  figure  6  for  n  =  0.707.  The  agreement  is  so 
close  that  the  two  sets  of  points  seem  to  define  a  single  curve,  the 
uppermost  curve  of  figure  6.  (Other  calculations  for  the  point 
defined  by  cp  =  60°,  plotted  in  fig.  7,  confirm  the  expectation 
that  this  close  agreement  becomes  progressively  impaired  as  n  is 
chosen  smaller  and  smaller.  The  error  is  within  5  percent  down  to 
n  =  0.32 .) 

Equation  (27)  is  the  correct  solution  for  the  "symmetrical 
wake  correction"  of  reference  1.  This  result  has  been  incorporated 
in  an  erratum  sheet  thereto. 


Flow  II 
(Special) 


Flow  II  has  more  academic  than  practical  significance.  It 
could  be  used  in  the  determination  of  the  loading  on  special  ailerons 
to  cancel  (approximately)  the  loading  in  the  wake  of  a  sweptback  wing 
with  the  ailerons  deflected.  These  ailerons  would  be  full  span,  or 
located  inboard  if  part  span,  with  a  vertical  fence  at  their  juncture 
to  isolate  the  two  ahead  of  the  region  of  trailing -edge  disturbance. 

Comparison  of  the  specifications  of  flow  II  and  flow  I  in 
figure  2  shows  that  II  is  antisymmetric  whereas  I  is  symmetric.  It 
is  therefore  necessary  to  use  antisymmetric  source -sheet  sectors  in 
the  superposition  process,  rather  than  the  symmetric  sectors  used 
for  flow  I.  The  change  is  effected  by  simply  changing  the  sign  of 

the  second  cosh-1  term  in  equations  (18)  and  (19) .  The  equation 
corresponding  to  equation  (20)  then  becomes 


2 

rtPi^l-O2 


5w 

aa. 


o12da1 


Op2-^ 


du 

36  = 


(28) 
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(The  Cauchy  principal  part  of  the  integral  is  to  he  taken.)  The 
integral  equation  that  results  on  setting  ba/bo  =  0  is  thus 


a 


a12-o2 


=  o  |a|  <  n 


(29) 


in  place  of  equation  (21) .  Equation  (29)  is  of  the  same  form  as 
equation  (21):  dw/do^  is  replaced  hy  O^bw/bO^  as  the  unknown 

function  of  C According  to  equation  (24),  a  solution  of 
equation  (29)  therefore  is 


3oEj*al  =  C  esc  0  (30) 

in  terms  of  the  function  0  defined  in  equations  (22). 

For  the  determination  of  u,  the  solution  equation  (30)  is 

2 

substituted  hack  into  equation  (28)  with  0^2  =  2—  (1  -  cos  0) 

2 

as  in  equation  (22).  There  is  obtained 


n 


du 

3a  “ 


C 


JC 


d0 


1  +  cos  0 


n  c  |a|  <  1 


This  expression  for  ba/bo  differs  from  equation  (25)  only  in 
the  lack  of  a  factor  0  outside  the  integral  sign.  Comparison 
with  equation  (26)  therefore  shows  that  the  result  of  the  inte 
gration  must  he 


du 

35  =  ‘ 


Cn‘ 


26  I  o  |  ^/(l-o2)(a2-n2) 


n  <  |0|  <  1 


(31) 
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This  function  is  readily  integrated  with  the  aid  of  the  elliptic - 
function  substitutions  of  appendix  C.  (The  absolute  value  sign  is 
temporarily  ignored.)  The  integral  from  c  to  1  can  be  expressed 
in  the  form 


u(o)  .  SB  °  cos-1  IUnV -2n2 

(uV  -  - 

where  the  sign  factor  <j/|(7|  has  been  appended  in  order  to  provide 
the  antisymmetry  of  u  in  ±0  indicated  by  equation  (31) .  The 
constant  C  is  established  by  the  condition  u(n)  =  un,  so  that 
finally  u 


u(a)  =  1^2  _2_  cos”1  .(1+n2 )q2 ~2n2  n  <  |a|  <  i  (32) 
*  lql  (i-n2)a2 


Equation  (32)  is  the  solution  for  the  u  velocity  in  the  w  =  0 
sectors  of  flow  II  (fig.  2).  As  before,  the  solution  is  confirmed 
by  comparison  with  the  solution  (equation  (12))  for  the  related 
flow  II'  of  figure  2,  with  attention  to  the  singularities  in  du/do. 
The  two  solutions  are  plotted  together  in  figure  6  for  n  =  0.707. 
Again,  the  agreement  is  such  that  the  two  sets  of  points  seem  to 
define  a  single  curve,  the  second  from  the  top  in  the  figure. 


Flow  III 

The  boundary  conditions  (fig.  2)  show  that  flow  III  cannot 
be  conical  like  flows  I  and  II,  but  must  be  quasi -coni cal.  That 
is,  the  velocity  components  are  not  of  the  form  u  =  F(y/x)  and 
w  =  G(y/x),  but  are  rather  of  the  form  u  =  xF(y/x)  and 
w  =  xG(y/x).  The  flow  is,  however,  antisymmetric  like  flow  II. 

A  flow  of  this  type  can  be  built  up  (compare  treatment  of  flows  I 
and  II)  by  superposing  elementary  source -sheet  sectors  that  have  a 
strength  variation  proportional  to  y  or,  what  amounts  to  the 
same  thing,  to  x<7. 

An  example  of  how  an  arbitrary  antisymmetric  quasi-conical 
w  velocity  distribution  of  the  form 


WC^)  =  XOpffcp) 


(33) 
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can  be  built  up  by  such  sheets  is  shewn  in  figure  8.  (Compare 
fig.  5,  section  A-A.)  The  strength  of  an  elementary  sheet  is 
characterized  by  the  strength  parameter 


The  incremental  u  velocity  induced  by  an  elementary  sheet  for 
which 


w  =  (constant)  Pyj_ 

=  (Af) 

can  be  derived  by,  standard  methods  (for  example,  reference  5). 
The  expression  is 


Au  =  u(o,a1) 


df 

do 


p —  u(a,o1)  dO-L 


(34) 


where 


U(a,a1) s 


al2 


Z  2x3/2 
*(1-0^) 


1  l+oo  ,  1-0  0 

(1+0-.0)  cosh  : - =— T  -  (1-0,  o)  cosh  - - i 

ai+CT  “ 


a1-a 


The  total  value  of  u  contributed  by  the  superposition  of  all  the 
contemplated  sheets  is  given  by 

nn 


u(x,o) 


55T  doi 


(35) 


Now  u(x,o)  is  specified  over  the  shaded  area  (|o|  <  n)  of 
flow  III  in  figure  2.  Thus,  equation  (35)  represents  an  integral 
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equation  for  the  unknown  function  df  The  kernel 

U(a>aq)  of  this  integral  equation  is  quite  complex.  The  second 
derivative  of  11(0,0^)  with  respect  to  0  has,  however,  the 
much  simpler  form 


^(0,0!)  _  40ax3 

do2  Jt  (a^-O2)2^!-^ 


This  simple  function  would  then  he  the  kernel  of  the  integral 
equation  obtained  by  differentiating  equation  (35),  provided  that 
the  right-hand  side  of  equation  (35)  could  be  differentiated  under 
the  integral  sign  without  regard  for  the  singularity  at  0q  =  a. 

It  is  clear  that  such  a  procedure  would  be  invalid  in  this  case 
because  for  |0|  <  n  the  kernel  would  have  a  second -order  pole 
within  the  range  of  integration. 


H.  Mirels  of  the  Lewis  laboratory  has,  however,  pointed  out 
that  if  the  differentiation  were  conducted  properly  the  result 
might  be  written  in  the  form 


=  i.  _££ _ 4a<^- _ d0T  (36) 

3  J0 

r 

The  symbol  signifies  that  the  indefinite  integral  evaluated 

JO 

at  a-]_=  0  is  to  be  subtracted  from  the  indefinite  integral 

evaluated  at  =  n.  Defined' in  this  way,  the  integration  can 

be  carried  out  without  difficulty,  inasmuch  as  the  integrand  is 
regular  in  the  neighborhood  of  the  limits. 

The  formulation  of  equation  (36)  is  a  special  case  of  a 
theorem  that  may  be  briefly  stated  as  follows :  If  the  function 
f(x,U  has  an  "integrable  singularity”  within  the  interval  of 
integration,  then,  under  rather  general  circumstances,  the  equation 


4 
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Pb 


P  b 


f(x,|)d| 


f(x,|)d| 

Sxn 


(D4) 


is  correct.  A  general  discussion  and  proof  of  this  theorem  is 
provided  in  appendix  D  by  F.  K.  Moore  of  the  Levis  laboratory. 

The  prescribed  condition  on  u  in  the  range  |c|  <  n  is 
(flow  III,  fig.  1): 


or 


u 


Ky 


Kxq 

P 


|C|  <  n 


Application  of  this  condition  to  equations  (35)  and  (36)  yields 
the  two  alternate  integral  equations  for  df  (q  ^ ) /dq^ 


Kq  = 


r  n 


A 


0(0,0! )do 


r  l°l  <  n 


0  = 


df  gl5agl 
dai  (o^-op 


(37) 


(38) 


Attention  will  be'  centered  on  solving  the  simpler  of  the  two 
equations,  (38). 
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It  Is  convenient  to  write  equation  (38)  in  the  form 


where 


P 

1  g(Oi)g1d01 

Jo 


g(o2) 


(39) 


(40) 


« 


is  regarded  as  the  unknown  function.  Then  equation  (39)  will  be 
solved  if  any  choice  of  gfa^)  can  be  found  that  will  cause  the 
integral  I  to  vanish. 


The  transformation 


-\ 


g(Oq) 

„  2 


1  " 


a  = 


Cn2h(8) 

n2  sin2  8 

n2q2 
q2-l 


J 


(41) 


converts  equation  (39)  to 


I 


C(l-q2)2 


P*/2 

i.  h(0)  sin  0  cos  9  d 0 

T  (sin2  8  +  q2  cos2  8)2 

JO 


(42) 


Certain  similarities  to  integrals  14  and  15  of  reference  11,  table  48, 
suggest 


4 


h(8)  =  tan  8 
h(8)  =  cot  6 


(43a) 

(43b) 


1  xi  n 
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as  possible  solutions  of  equation  (38)  or  (39) .  Both  functions 
are  found  upon  substitution  and  integration  to  satisfy  the 
equation.  Any  linear  combination  of  solutions  (43)  is  also 
clearly  a  solution.  Thus,  the  integral  equation  (38)  does  not 
have  a  unique  solution. 

The  later  identification  of  equation  (43a)  as  the  physically 
correct  solution  is  anticipated  at  this  point.  Then  the  u  velocity 
in  the  w  =  0  regions  of  flow  III  may  be  evaluated  as  follows. 
Equation  (36)  can  be  written 


(44) 


where  I  is  the  integral  defined  in  equation  (39)  or  equation  (42). 
With  the  solution  equation  (43a)  substituted  in  equation  (42),  the 
integral  I  has  the  value  (reference  11,  table  48,  integral  14) 


I  =  C(l-q2)2  JL  n  < |a| <  1  (45) 

4q 

But  according  to  equations  (41) 


whence 


I 


Cn4 

o(cr2-n2)5/2 


n  <  |a|  <  1 


and,  upon  returning  to  equation  (44), 


fu 


4 

Cnx _ 

P(a2-n2)3/2il^ 


n<  f a[  <  1 


(46) 
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The  u  velocity  in  the  region  n  <_o  <  1  results  upon  inte¬ 
grating  equation  (46)  twice  between  the  limits  a  and  1.  This 
integration  may  he  effected  with  the  aid  of  the  elliptic -integral 
substitutions  of  appendix  C,  and  the  result  is 


The  value  of  C  is  established  by  the  condition  that 
u(n,x)  =  Ky  =  Knx/6.  The  negative  range  of  O  (-1  <_a  <_  -n)  is 
taken  care  of  by  introduction  of  the  antisymmetry  factor  o/\o\. 
These  operations  yield  the  final  result 


n  <  |a|  <  1 


(47) 

Equation  (47)  is  the  solution  for  the  u  velocity  in  the  w  =  0 
sectors  of  flow  III  (fig.  2).  This  equation  results  from  the  choice 
of  equation  (43a)  among  the  solutions  of  the  integral  equation  (39) . 
The  correctness  of  this  choice  is  determined  by  comparison  of 
equation  (47)  with  the  solution  equation  (15)  of  the  related  flow  III 
of  figure  2.  Both  equations  exhibit  a  half -order  singularity  in 
du /bo  at  O  =  n+  together  with  du/da  =0  at  a  =  1. 

(Equation  (43b),  on  the  other  hand,  leads  to  a  different  behavior.) 
Equations  (47)  and  (15)  are  plotted  together  in  figure  6  for 
n  =  0.707.  Again,  as  for  flows  I  and  I'  and  II  and  II'  the  agreement 
is  such  that  the  two  sets  of  points  lie  on  a  single  curve,  the 
bottom  curve  in  the  figure.  The  degree  of  agreement  as  a  function 
of  the  parameter  n  for  the  ray  defined  by  cp  =  60°  is  shown  in  fig¬ 
ure  7.  The  error  is  within  5  percent  down  to  n  =  0.4. 


U«T,X)  .  j§:  f 


_  _  2 

\a\  E(cp,k)-n2F(cp,k)  -  sin  2Cp 

Efe,k\  -n2F(5-,k 
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Flow  IF 

The  boundary  conditions  (fig.  1)  show  that  this  flow  is  quasi- 
conical  like  flow  III.  That  is,  the  velocity  components  are  of  the 
form,  u  =  xF(y/x)  and  w  -  xG(y/x).  The  flew  is,  however,  symmetric 
like  flow  I.  A  flow  of  this  type  can  be  built  up  (compare,  treat¬ 
ment  of  flows  I  and  III)  by  superposing  elementary  source -sheet 
sectors  that  have  a  strength  variation  proportional  to  x.  These 
sheets  are  arranged  in  the  same  manner  as  the  corresponding  constant - 
strength  sheets  employed  for  flow  I.  (See  section  A-A,  fig.  5.) 

Let  the  source  strength  of  an  elementary  sheet  sector  be 
specified  by 

w  =  (Af)  x 

where  Af  is  a  constant  for  a  particular  sheet  and  the  sector 
extends  from  a  =  -  Op  to  o  =Op  The  incremental  u  velocity 

induced  by  such  a  sheet  has  been  calculated  by  the  method  of 
reference  5  and  the  result  is 
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(The  Integral  was  not  evaluated  inasmuch  as  only  the  x -derivatives 
of  W  will  he  needed  in  the  applications  herein.) 

The  individual  sheet -strength  parameter  Af  can  be  expressed 
as 


Af 


dc^ 


so  that 

df 

do^ 

Au  =  -  -jf  W (x,©,^)  d^ 


The  total  value  of  u  contributed  by  the  superposition  of  al 1  the 
contemplated  sheets  is  then 


n 

/  k  1  df  , 

u(x,a)  =  -  p  (x,^)  d a±  (49) 

Jo 

The  development  beyond  thiB  point  runs  similarly  to  that  for 
flow  IU.  Inasmuch  as  u(x,o)  is  specified  over  the  shaded  area 
(|Ol<  n)  of  flow  IV  of  figure  2,  equation  (49)  is  an  integral 
equation  for  the  unknown  function  df(a1)/dO]L.  The  kernel 

W(x,a,a3_)  is  disagreeably  complex,  but  its  second  derivative 
with  respect  to  x. 


d2 

^w(x,a,ffL) 


y 


-4CT4a-L 

jtx^-c?)2  ]/i^ 


is  relatively  simple.  The  theorem,  equation  (D4),  for  differenti¬ 
ation  under  the  integral  sign  is  therefore  employed  to  give 


df  _ 

d°l  nxiOj2-#)2 


(50) 
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rn 


The  integration  symbol  ^ 

Jo 

discussed  under  flow  III. 


has  the  special  significance  previously 


The  prescribed  condition  on  u  in  the  range  |a|  <  n  is 
(flow  IV,  fig.  2)  “ 


or 


>1 

u  =  kx 


Application  of  this  condition  to  equation  (50)  yields  the  integral 
equation  for  df (0^)/dC^ 


Jo 


|C|  <  n 


(51) 


This  equation  may  be  compared  with  the  corresponding  equation 
for  flow  III  (equation  (38)).  Both  equations  are  of  the  form  of 
equation  (39),  but  with  different  expressions  for  the  unknown 
g(a1).  Thus  the  method  of  solution  for  equation  (39)  will  apply 

for  equation  (51),  with 


«<V  -  if 


(52) 


Once  again  the  proper  solution  for  g(o^)  =  Cn2h(0)  is 

found  (by  later  comparison  of  flows  IV"  and  IV')  to  be  given  by 
equation  (43a).  This  solution  is  now  employed  in  equation  (50) 

in  the  range  n  <  |0|  <1  for  which  c^u/cbc2  does  not  vanish. 
The  integral  has  the  value  given  in  equation  (45),  and  equation 
(50)  becomes 
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CnV 


^u 

ax2  L  .  "  ,.3/2 


y  Px(o2-n2)  ^jl-o2 


n<  |a|  <1 


(53) 


Now  u  is  known  to  be  of  the  foiro  (quasi -conical  flow) 


so  that 


u  =  Pyf(c) 


"\ 


=  -  cfif'io)  =  g (a),  say 


(54) 


J 


Thus  "by  comparison  of  equation  (53)  with  the  last  of  equations  (54) 


g’(CT)  = 


CnV 


P  9  3/2 
P(a2-n2) 


fl. 


The  elliptic -function  substitutions  of  appendix  C  are  convenient 
for  performing  the  integration  for  g(o) ,  with  the  limits  a 
and  1.  Then  f'(a)  =  -g(a)/i2  may  in  turn  he  integrated  to 
yield  u.  The  result  is 


u(c,x)  = 


Cn4x 


P(l-n2) 


[> 


F(cp,k)-E(co,k) 


where 


cp  =  sin”'1'/^ 

P 

k  =  \jl- n2 
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The  constant  C  is  established  by  the  condition  that 

u(n,x)  =  kx 


so  that  finally 


u(0,x)  =  K  x 


F(cp,k)  -  E(cp,k) 


(55) 


Equation  (55)  is  the  solution  for  the  u  velocity  in  the 
w  *  0  sectors  of  flow  IV  of  figure  2.  This  equation  results 
from  the  choice  of  equation  (43a)  among  the  solutions  of  the 
integral  equation  (39) .  The  correctness  of  this  choice  is  deter¬ 
mined  by  comparison  of  equation  (55)  with  the  solution  equation  (17) 
of  the  related  flow  IV*  of  figure  2.  Both  equations  exhibit  a 
half -order  singularity  in  du/da  at  0  =  n+  together  with 
du/da  =0  at  0  =  1.  (Equation  (43b),  on  the  other  hand,  leads 
to  a  different  behavior.)  Equations  (55)  and  (17)  are  plotted 
together  in  figure  6  for  n  =  0*707.  Once  again,  as  for  the  pairs 
of  flows  I  and  1 1 ,  II  and  II  * ,  and  III  and  III 1 ,  the  agreement  is 
sufficiently  close  so  that  the  two  sets  of  points  seem  to  define  a 
single  curve,  which  is  the  curve  next  to  the  bottom  in  the  figure. 
The  degree  of  agreement  as  a  function  of  the  parameter  n,  for  the 
point  defined  by  cp  =  60°,  is  shown  in  figure  7.  The  error  is 
within  5  percent  down  to  n  =  0.27. 


DISCUSSION 

In  the  foregoing  derivations  principal  attention  has  been 
devoted  to  the  solution  for  the  u  velocity  (or  the  pressure  dis¬ 
turbance)  in  the  regions  in  which  it  is  unknown.  The  solution  for 
the  w  velocity  (or  the  upwash)  in  the  regions  in  which  it  is 
unknown  has  also  been  obtained,  although  for  brevity  it  has  been 
omitted  herein  (except  in  the  single  example,  appendix  B).  Com¬ 
parisons  of  the  singularities  in  the  w  velocity  between  flows  I 
and  If,  III  and  III1,  and  IV  and  IV1,  respectively,  were  made  in 
addition  to  the  u  velocity  comparisons  mentioned  earlier.  These 
additional  comparisons  again  showed  very  close  agreement  between 
each  flow  and  its  related  flow,  in  support  of  the  correctness  of 
the  choice  in  each  case  among  the  solutions  of  the  integral 
equations  for  flows  I,  IH,  and  IV. 
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A  primary  aim  of  this  report  has  been  to  provide  the  analytical 
basis  for  approximating  in  simple  fashion  the  effect  of  subsonic 
trailing  edges  on  damping  in  pitch  and  roll  for  thin  sweptback  wings 
in  a  supersonic  stream.  (See  reference  3  for  detailed  application.) 
Accordingly,  the  main  emphasis  has  been  placed  on  the  "desired  flows" 

I  to  IV.  The  more  easily  solved  "related  flows"  I'  to  IV'  were  found, 
however,  to  approximate  flows  I  to  IV  in  the  regions  of  interest 
much  better  than  had  been  anticipated.  (See  figs.  6  and  7.)  It  is 
clear  from  the  comparison  that  a  moderate  modification  of  the  boundary 
condition  near  the  left-hand  Mach  line  (fig.  2)  -  the  distinguishing 
feature  between  the  unprimed  and  primed  flows  -  will  scarcely  affect 
the . flow  in  the  region  near  the  right-hand  Mach  line.  Thus  such  a 
modification  may  freely  be  made  to  simplify  a  given  problem,  and  the 
solution  of  the  modified  flow  will  apply  to  the  desired  flow  with  " 
engineering  accuracy.  In  particular,  the  related  flows  herein  may 
be  used  in  place  of  the  desired  flows  if  trigonometric  functions  are 
preferred  to  elliptic  integrals. 

A  general  method  for  the  calculation  of  flows  (fig.  9)  of 
which  flows  I'  and  IV'  are  special  cases  is  included  in  a  paper  by 
Goodman  (reference  12)  after  the  bulk  of  the  present  work  was 
completed.  The  method  and  its  applications  have  been  simplified  * 

and  extended  by  Mirels  (reference  13).  The  general  solution  for 
the  unknown  u  velocity  in  the  flows  specified  in  figure  9  is  given 
in  elegant  form  in  reference  13.  This  result,  originally  obtained 
by  means  of  a  doublet  distribution,  may  be  derived  by  means  of  a 
development  of  the  present  source -distribution  method.  The  unknown 
w  velocity  may  be  obtained  as  well.  The  details  are  given  in 
appendix  E  as  a  matter  of  interest. 


CONCLUDING  REMARKS 

The  flow  over  a  sweptback  wing  may  be  obtained  by  superposing 
on  a  basic  delta-wing  flow  additional  flows  to  cancel  the  lift 
outside  the  boundaries  defining  the  sweptback  wing.  The  cancel¬ 
lation  flow  for  the  trailing  edge  modifies  the  lift  in  a  region 
ahead  of  the  trailing  edge  if  the  component  of  the  stream  velocity 
normal  to  the  edge  is  subsonic.  For  angle  of  attack  the  principal 
part  of  this  cancellation  flow  is  conical;  for  rolling  and  for 
pitching  the  principal  part  is  quasi -conical .  The  derivation  of 
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these  and  several  related  conical  and  quasi -conical  flows  has 
heen  carried  out  in  this  report.  For  each  case  the  problem  was 
formulated  as  an  integral  equation.  Some  of  the  results  are 
applied  to  damping  in  roll  and  in  pitch  in  another  report. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  February  3,  1950. 
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APPENDIX  A 
SYMBOLS 

The  following  symbols  are  used  in  this  report: 

C  boundary  for  line  integral,  as  specified  in  text;  also 

undetermined  constant 

Cp  pressure  coefficient  (proportional  to  u) 

E(cp,k)  incomplete  elliptic  integral  of  second  kind  with  amplitude 
cp  and  modulus  k 

F(cp,k)  incomplete  elliptic  integral  of  first  kind  with  amplitude 
<P  and  modulus  k 

f(0)  function  defined  in  equation  (22) 

f (a)  function  defined  in  equation  (54) 

f(a^)  function  defined  in  equation  (33) 

g (a)  function  defined  in  equation  (54) 

gfcl)  function  defined  in  equation  (40)  for  flow  III,  in 

equation  (52)  for  flow  TV 

h(@)  function  defined  in  equation  (41) 

I  integral  defined  in  equation  (39) 

J  integral  defined  in  equation  (Bl) 


M  Mach  number,  ratio  of  stream  velocity  to  velocity  of 

sound  in  free  stream 

n  value  of  o  along  side  edges  of  shaded  triangular  regions 

in  figure  2 


S  area  for  surface  integral,  as  specified  in  text 
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11(0,^) 

u,v,w. 

u0 

Au 

Y 

x,7,z 

a 

P 

5 

9 

K 

i>*\ 

a  = 

t  = 

cp  = 

I, II 
III,IY 

m', iv* 


function  defined  following  equation  (34) 

disturbance  velocity  components  along  x-,  y-,  and 
z-axes,  respectively 

constant  value  of  u  over  prescribed  area 
increment  in  u 


free -stream  velocity 

function  defined  following  equation  (48) 

Cartesian  coordinates:  x-axis  parallel  to  free -stream 
direction;  y-axis  horizontal  and  toward  right,  looking 
upstream;  z-axis  vertically  upward 

angle  of  attack 

Vm2-i 

defined  in  equation  (22) 

defined  in  equation  (22)  for  flows  I  and  II,  in 
equation  (41)  for  flows  III  and  IY 


constant  of  proportionality  (fig.  2) 


oblique  coordinates  measured  parallel  to  downstream  Mach 
lines,  defined  by  equations  (2) 

Py/x 


(See  fig.  3  for  geometric  significance.) 


desired  flows  designated  in  figure  2 


related  flows  designated  in  figure  2 
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r* 


special  integration  sign  defined  following  equation  (36) 


A  primed  function  signifies  the  first  derivative  of  the 
function  with  respect  to  the  independent  variable. 
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APPENDIX  B 


EVALUATION  OF  UPWASH  IN  REGION  B  OF  FLOW  I' 


Note  that  just  Before  equation  (4)  there  is  the  definition 


Now  J  is  a  known  function,  which  was  evaluated  in  the  course  of 
obtaining  equation  (10) .  Thus  equation  (Bl)  is  an  Ahel  integral 
equation  for  the  unknown  function  dw/2 be,.  Its  solution  (refer¬ 
ence  6)  Is 


dv 

^1 


=  M 


a 


VT 


(B2) 


With  the  value  of  J  for  flow  I*  (equation  (6))  this  is,  provided 
that  the  integral  on  the  right-hand  side  is  continuous  (reference  7), 

VE=i 

riiA 


dw  Muo  a 


Mu0 

*  ^idi-niA) 

Upon  dropping  the  subscripts  and  converting  to  Cartesian  coordinates 
by  means  of  equations  (2), 

Bw  =  Pu0,l  2(l+n) 

&  Jt  ^(x+PyKnx-Py) 


(B3) 
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The  indefinite  integral  is 


w 


Pu0  |2(l+n) 
n  \j  n 


cosh"1  2ni-Py(l-n) 
Py  (1+n) 


Pu0  /2  (1+n) 

n  V  n 


cosh 


-1 


2n-q(l-n) 
|a|  (1+n) 


(B4) 


Equation  (B4)  gives  v  =  0  along  the  |-axis  (O  =  -1)  and  along 
the  line  q  =  (o  =  n),  as  it  should  for  flow  I’.  The  constant 
of  integration  is  therefore  zero  and  equation  (B4)  is  the  required 
solution  for  the  upwash  velocity  w. 
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APPENDIX  D 

DIFFERENTIATION  OF  INTEGRAIS  OF  SINGULAR  FUNCTIONS 
By  Franklin  K.  Moore 

The  problem  arising  in  the  differentiation  of  equation  (35) 
may  he  stated  in  the  following  general  terms :  It  is  required  to 
differentiate,  with  respect  to  a  parameter  of  the  integrand,  the 
definite  integral  of  a  function  that,  though  integrable,  has  a 
singularity  inside  the  range  of  integration  such  that  the  definite 
integral  of  the  derivative  is  not  convergent .  It  is  shown  herein 
that,  subject  to  certain  restrictions,  this  differentiation  can  be 
carried  out  in  a  simple  manner,  avoiding  consideration  of  the 
singularity. 

The  integral  to  be  differentiated  can  be  written 

'"’b 

I(x)  =  f(x,!)d* 

Ja 


It  will  be  supposed  that  f  (x, |)  has  a  singularity  at  some  point 
a  <  i  <  b  and  is  a  function  such  that 


(l)  An  indefinite  integral  F(x, I) 


found . 


f 

[  f(x,t)d£  may  he 


(2)  The  function  l(x)  is  a  convergent  improper  integral  that 
can  be  written  l(x)  =  F(x,b)-F(x,a) .  (See  reference  16,  para¬ 
graph  169.)  Under  these  restrictions,  f(x)  may  have  an  "integrable 
singularity"  within  the  interval  of  integration.  For  example, 
f (x,t )  =  logg  | x—  1 1,  a  <  x  <  b,  meets  the  foregoing  requirements 
(1)  and  (2) . 


It  is  required  to  carry  out  the  following  differentiation: 
b 

^  \ ~ ^  _/  _\  S 


I’(x)  =  /  f(x,|)d£  =  F(x,b)  -  F(x,a) 


(Dl) 
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or 


i’M  = 


j 


r  a 


^•f(x,|)d|-  ^f(x,f)d| 


When  the  convention  introduced  in  equation  (36)  is  followed, 
equation  (D2 )  can  he  written 


(D2) 


P  - 

I’to-i  £f(x,tMl  (D3) 


The  result  (D3)  may  be  stated  as  follows:  Subject  to  the 
restrictions  (l)  to  (4)  on  the  behavior  of  the  integrand  near 
the  limits,  the  differentiation  with  respect  to  a  parameter  to 
the  integrand  of  a  convergent  improper  integral  may  be  accomplished 
by  formally  integrating  the  derivative  of  the  integrand,  as  though 
its  singularity  were  not  present. 

It  is  clear  that 

P 

I"M  =  |  ~-f(x,|)d£  (D4) 

Ja 


provided  that  requirement  (3)  is  replaced  by  the  restriction  that, 
near  |  =  a,b 


ctx' 


FMn 


a 


Equation  (D4)  is  the  relation  used  in  this  report. 

The  foregoing  proof  may  easily  be  extended  to  the  case  where 
a  and  b  are  functions  of  x,  subject  only  to  the  additional 
requirement  that 
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APPENDIX  E 


GOODMAN -WHEELS  GENERALIZED  CANCELLATION  FLOW 

The  generalized  cancellation  flow  introduced  in  reference  12 
and  extended  in  reference  13  is  specified  by  the  boundary  condi¬ 
tions  of  figures  9(a)  or  9(b),  depending  on  the  position  of  the 
line  i]  =  f(|).  Consider  first  figure  9(a).  The  flow  can  be 
represented,  as  before,  by  sources  distributed  along  the  surface 
in  proportion  to  w.  In  the  Mach  coordinates  of  reference  5,  the 
surface  potential  is 


cp(bn) 


\jn-n i 


'chi) 


vnr 


region  B 


where  (|,r))  is  located  at  P  in  region  B.  In  what  follows  <P 
will  play  essentially  the  role  that  u  played  in  the  solution 
of  flows  I'  to  IV' . 

Consider  next  the  case  of  figure  9(b).  By  virtue  of  the  left- 
hand  u  =  0  region  (reference  5),  the  region  of  integration  for 
cp(t,Ti)  will  not  include  the  area  left  unshaded  in  the  upper  corner. 
It  is  found  that  with  the  present  notation  the  limits  of  integration 
are  unchanged  from  the  case  of  figure  9(a).  Thus  equation  (El) 
applies  equally  well  for  the  case  of  figure  9(a)  or  9(b). 

Solution  for  u  in  region  A.  -  Define 


G(tJi  ) 


wd^ 

VhT 


so  that 


<?(!>*])  =  - 


region  B 
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Equation  (E3)  is  an  Abel  integral  equation  for  J.  Its  solution 
is  (references  6  and  7) 


^l-2 


Now  let  the  general  point  (|,t))  he  located  in  region  A 
at  Q.  The  only  change  in  the  expression  for  <T>  is  in  the  upper 
limit  for  ti^*  This  upper  limit  is  now  at  the  intersection  of 

the  line  PQ  with  the  line  r\  =  g(|).  Accordingly, 


Cp(l>Tl) 


)dq 


1'  ‘1 


The  value  of  J  already  obtained  in  equation  (E4)  applies  here 
as  well  as  earlier.  Upon  making  the  substitution, 


Vilrt)  =  Z 


"1g(D  Pi! 

- _<L  SLLils 

\jrip  z 

f  (I)  Jf(l) 


region  A 


_^Tl_ 


<m»,z.)d-z 


region  A  (E7) 


where  the  sign  |  designates  the  finite  part  of  the  integral 

according  to  Hadamard's  definition  (reference  17). 
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The  factor  l/  which  does  not  contain  z,  may  be 

included  in  a  common  integrand,  and  the  finite -part  sign  may  be 
removed  to  the  outside  of  the  double  integral: 


«p(!>n)  =  - 


1 

2rt 


"'g(l) 

,f(|) 


°nl 

gP(l,z)dz 

/\Tri^r  (ri1-z)3^ 

f(t)  V 


(E8) 


The  order  of  integration  may  now  be  reversed  with  the  aid  of  the 
scheme  for  the  limits  sketched  in  figure  10: 


¥(1,11)  =  - 


1_ 

2it 


nsd) 


f(i) 


cp(  |,z)dz 


e(S) 


dti- 


l-Tll  (ril“z)3/^ 


(E9) 


Evaluation  of  the  finite  part  of  the  inner  integral  yields  the 
final  result 


<p(  £>T|) 


Gm. 


Pg(|) 


^f(S) 


ffdjzHz 
d-z)  ^gd)-z 


region  A  (E10) 


Equation  (E10)  gives  the  surface  potential  in  region  A  that  results 
from  a  prescribed  distribution  of  surface  potential  in  region  B* 

Equation  (E10)  is  equivalent  to  equation  (15)  of  reference  15, 
which  applies  to  the  mirror  image  of  figure  S.  (The  same  equation 
was  obtained  for  a  restricted  situation  in  reference  18  (equation  (61d) 
therein).)  According  to  reference  13,  differentiation  yields 
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» 


u(t,n) 


..  ^SH) 


5t 


■g  (t) 


f(l) 


u( £,z)dz 


(t-z)  A^g(|)-z 


1-g'd) 


2rtP  ^ti-gd) 


ur 


*e(|) 


d) 


Pu  +  V 


dz 


region  A 

(Ell) 


which  is  equation  (17b)  therein.  This  equation  constitutes  a 
general  solution  for  the  disturbance  in  region  A  regardless  of 
whether  the  Kutta  condition  is  imposed  along  the  line  r|  =  g(|). 

It  is  shown  in  reference  13  that  imposition  of  the  Kutta  condition 
along  the  line  T]  =  g(|)  modifies  the  sidewash  distribution  v 
in  such  a  way  as  to  cause  the  second  integral  to  vanish.  (If  the 
boundary  tj  =  g(£)  of  fig.  9  slopes  everywhere  toward  the  left, 
instead  of  toward  the  right  as  shown,  v  is  determined  uniquely 
by  the  known  u  distribution,  and  the  Kutta  condition  may  not 
be  imposed.) 


The  Kutta  condition  version  of  equation  (Ell)  (that  is,  with 
the  second  integral  set  equal  to  zero)  may  be  obtained  more  directly 
for  case  (b)  of  figure  9.  In  this  case  (including  the  special 
case  f (I )  =  0) ,  u  may  be  written  in  place  of  cp  and  dw/cfac^ 

in  place  of  w  in  equations  (El)  through  (E10);  the  result  is 
given  in  the  new  equation  (E10) .  This  approach  at  first  parallels 
and  then  extends  that  of  equation  (3)  and  the  following  equations 
in  the  text. 


Solution  for  w  in  region  B.-  It  is  worthwhile  to  complete 
the  solution  of  the  Goodman -Mirels  cancellation  flow  (fig.  9)  by 
evaluating  the  upwash  velocity  w  in  region  B.  This  knowledge 
of  w  will  be  useful  in  determining  the  downwash  field  in  the 
general  vicinity  of  the  wing  -  in  particular,  within  a  chord  length 
behind  the  trailing  edge  -  by  means  of  the  cancellation  technique. 
(A  line  vortex  method  (references  19  and  20)  yields  the  downwash 
farther  back  more  simply,  and  simple  means  are  known  (cited  in 
references  19  and  20)  for  determing  the  downwash  immediately  behind 
the  trailing  edge . )  * 
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The  solution  is  started  "by  inverting  (solving)  the  Abel 
integral  equation  (E2)  with  the  aid  of  references  6  and  7: 


w 


G(i)l) 


J(!,r|i)a£ 


(E12a) 


J(l  ^T]1)d| 

(li-l)372 


(E12) 


The  solution  for  J  has  already  been  determined  in  equation  (E4) . 
Changing  z  to  q  therein  gives 


i  d 
*  ^ 


p\ 


Jf  (£) 


<P(^»h)dn 

Vv^ 


(E13a) 


_1 

2it 


^i 

f(D 


l>Ti)dri 

(Tll-tl)^ 


(E13) 


Substitution  of  equations  (E13a)  and  (E13)  into  equations  (E12a) 
and  (E12 ) ,  respectively,  yields 
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»(ti,m>  =  *  5 


M  5 


3f 


*L 


G(tii) 


4  t 


\ji i“l  l^i 


^l 


fit) 


Cp(l »  ti  )  Art 


M_ 

4it 


li 


(1 1-D 


372 


region  B 


(E14a) 


region  B 

(tj  -n)3/2 

f(|)  1 


(E14) 


By  comparison  with  equation  (10)  of  reference  13,  equation 
(E14)  may  he  interpreted  as  the  upwash  velocity  w  due  to  a 
given  distribution  of  doublets  (or  of  vorticity)  over  a  prescribed 
area.  The  area  of  integration  is  shown  cross-hatched  in  figure  11 
for  the  same  two  configurations  (a)  and  (b)  considered  in  figure  9. 
Note  that  the  regions  A'  and  B',  which  lie  in  the  zone  that  influences 
the  point  are  not  included.  This  result  implies  that  for 

the  point  the  integral  of  the  doublets  in  A’  exactly  cancels 

the  integral  of  the  doublets  in  B*.  The  situation  here  with  respect 
to  doublets  presents  an  interesting  parallel  to  Evvard's  original 
discovery  with  respect  to  sources  in  connection  with  his  membrane 
concept.  (See  reference  21.) 


Equation  (E14a)  is  an  alternate  form  of  equation  (E14)  that 
may  yield  simpler  integrations  in  some  cases.  Note  that  the  | 
(finite  part)  operation  is  avoided.  Similar  considerations  apply 
to  the  alternate  forms  (E12a)  and  (E12),  (E13a)  and  (E13). 

A  particularly  simple  result  for  dw/dx  can  be  obtained 
when  the  Kutta  condition  is  imposed.  Thus,  apply  to  equation  (E14) 
a  procedure  analogous  to  that  by  which  dcp/Bx  is  obtained  from 
cp  in  reference  5.  The  result  is 
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The  first  line  integral  vanishes  because  cp  =  0  along  ah,  and 
the  second  line  integral  may  be  expressed  differently: 
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(E15) 


where  |2  =  g(Ti  )  • 
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The  line  integral  may  he  put  in  more  convenient  form  by  means 
of  an  integration  by  parts.  The  final  result  may  be  written 
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(E16) 


The  line  integral  is  the  same  as  the  one  in  equation  (Ell).  Here 
again,  then,  the  line  integral  vanishes  if  the  Kutta  condition  is 
imposed  along  the  line  \  =  G(r)') . 
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Figure  1.  -  Sweptback  wing  showing  basic  delta-wing  load 

distribution  along  section  A-A  and  approximate  cancellation 
of  this  load  behind  wing  by  superposition  of  special  flows. 
Special  flows  are  plotted  with  reversed  sign.  Cross- 
hatched  areas  represent  induced  changes  in  loading  ahead 
of  trailing  edge. 
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Figure  3.  -  Data  for  flow  I'# 


Figure  4.  -  Data  for  flow  II* • 
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Figure  5*  -  Superposition  of  elementary  source-sheet 
sectors  to  build  up  flow  I.  (Schematic) 
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Figure  8.  —  Superposition  of 
sectors  to  build  up  flow 
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